Realizing the driven nonlinear Schrodinger equation with stationary light 
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A tunable optical waveguide system to simulate the nonlinear Schrodinger equation (NLSE) with 
a periodic lattice potential in an out-of-equilibrium scenario is proposed. The system is formally 
similar to other nonlinear systems such as Bose-Einstein condensates in optical lattices or electro- 
magnetic fields propagating in conventional Kerr-nonlinear optical media. We discuss the advantages 
of the proposed system compared to these systems and investigate the transport properties in both 
the linear and nonlinear regimes. Multistable behaviour and a novel horizontal peak in the trans- 
mission spectrum is observed in the nonlinear regime, the latter near a band-gap region introduced 
by the periodicity. The proposed model provides an excellent testbed to study out-of-equilibrium 
NLSE dynamics due to the ease of introducing a coherent source (an input laser field), the extra 
control over the sign and strength of nonlinearity and the freedom in choosing the shape and depth 
of the underlying potential. As we discuss in detail, our proposal is within reach of experimental 
implementation. 

PACS numbers: 42.50.Gy,05.60.Cd,42.65.Pc 



I. INTRODUCTION 

Nonlinear Schrodinger equation (NLSE) arises natu- 
rally in many areas of physics. For example, it describes 
propagation of an electromagnetic field in a nonlinear 
medium [lj] or the interacting ground state of a Bose- 
Einstein condensate where the equation goes by the name 
of Gross-Pitaevskii equation [2[. Recently, it has been 
shown that the equation also governs the dynamics of 
stationary photons in a nonlinear optical fibre setup 
where strong interaction between photons are induced 
using enhanced Kerr nonlinearities obtained by electro- 
magnetically induced transparency (EIT) [4[. 

One of the interesting properties to study in systems 
governed by NLSE is nonlinear transport properties. Ex- 
amples include nonlinear transport of Bose-Einstein con- 
densates and soliton propagation in nonlinear op- 
tics p], Q . Photonic transfer in a nonlinear optical fibre 
set-up based on EIT has also been studied [10], but un- 
like in the previous examples, the study focuses on the 
regime of strongly interacting few-quanta without a po- 
tential. This nonlinear optical fibre system is an example 
of driven-dissipative systems that has attracted a lot of 
interest recently (see, for example, [lll - [l3| ) , and is natu- 
rally operated out of equilibrium. It is therefore an ex- 
cellent system to study the nonequilibrium physics such 
as the transport properties. 

In the present work, we study, in the semi-classical 
limit, the transport of photons inside a strongly nonlin- 
ear waveguide of length L depicted in Fig. [1] (a). The 
semi-classical regime is likely to be the starting point of 
experimental investigations which could be used not only 
to study the physics in the regime, but also to extract rel- 
evant parameters of the system for example. 
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Figure 1: a) Schematic diagram of the system, describing the 
quantum transport in an effective polaritonic lattice inside a 
waveguide, b) A Level diagram showing various detunings 
and coupling between the fields and levels. 



Taking advantage of the tunability of quantum optical 
systems based on EIT, one can vary the effective mass 
and the interaction parameters over a wide range of val- 
ues. In particular, we introduce a periodic lattice po- 
tential in the NLSE following the proposal of [14]. The 
proposed system thus allows the study of transport prop- 
erties in a specific periodic nonlinear system which has 
not been studied in detail before. As an example of a non- 
linear phenomena we investigate bistability behaviour in 
the transmission spectrum that have been well-studied in 
other nonlinear systems. 

The system studied in this work is formally similar 
to a Bose-Einstein condensate (BEC) in an optical lat- 
tice Il5| and also systems studied in nonlinear optics 
PI B Ell • It therefore provides a new versatile platform to 
study transport properties of the NLSE equation arising 
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in different physical systems, with several advantages. To 
name a few before we discuss them in more detail later, 
the platfom allows more efficient measurement of the out- 
put states optically, a more natural implementation of a 
coherent source (the drive) and extra control in the shape 
of the potential. 

Finally, we note that the nonlinear optical fiber sys- 
tem allows the study of quantum transport properties 
of interacting bosons 10]. We plan to investigate such 
quantum transport properties with a lattice potential in 
future work. 



II. SYSTEM 

The nonlinear waveguide system investigated in this 
work is based on the four level EIT scheme. A schematic 
diagram of the four level configuration and the coupling 
between the fields and levels are shown in Fig. [TJb). The 
quantum fields E± couple the ground state \a) and an 
excited state \b). A metastable state |c) is coupled to \b) 
through two classical counter-propagating fields of Rabi 
frequencies Q±. The nonlinearity is introduced in the sys- 
tem by adding an additional state \d), which is coupled 
to |c) throu gh t he qua ntum fields E±. Earlier related 
proposals [^Jiol [13 . fl7L fl8| were based on a hollow core 
fiber filled with ultracold atoms [l9l-[24). However, an- 
other ID waveguide system, tapered fiber, coupled with 
an ensemble of atoms has been demonstrated recently 
[25l [26[ and the authors feel that the both are promising 
systems to realize the particular Kerr-nonlinear photonic 
system considered in this work. 

Following the methods described in 0, [To|, we define 
dark-state polariton operators which, in the slow light 
limit, are written as = gy/2wnoE±/Q±. Here, no is 
the mean density of an ensembles of atoms coupled to 
the waveguide and g is the atom-field coupling strength. 
For simplicity, we restrict ourselves to the case where the 
Rabi frequencies of the control fields are equal: Q± = Q. 

Introducing the symmetric and anti-symmetric combi- 
nations of the polariton operators \I/ = and A = 



^■7= — , it is easily seen that, in the limit of large optical 
depth, A adiabatically follows ^: A = -ik L coh ^d z ^ . 
Here, m is the complex effective mass, htlr is the real 
part of m, ko is the wave vector of the lattice potential 
and L co h is the coherence length introduced in [10]. The 
specific forms of these parameters, and other effective in- 
teraction parameters to be introduced below, in terms of 
bare optical parameters are provided in Appendix [Aj 

As shown in [14], we introduce a polaritonic lattice 
potential in the above setup by creating a periodic mod- 
ulation in the atomic density in the ground state \a): 
n a — n o + n i cos 2 (k^z). The dynamics of the sym- 
metric polariton operator is governed by the nonlinear 
Schrodinger equation with a lattice potential, 

idt® = — —d*V + 2Vcos 2 (k z)^ + 2x^ 2 , (1) 
2m 



where, V is an effective lattice depth and x is an effective 
polariton-polariton interaction strength. After substitut- 
ing in ^f(z,t) = i/j(z, e)e~ %et 1 the above equation reduces 
to a dimensionless form: 



-d\i) + — [V - e + V cos(2z) + 2x^ty] ^ = 0, (2) 

where, z, ip and V has been scaled by 1/fco? V^o and the 
recoil energy Er = /(2m^), respectively. In terms of 
the right and left propagating fields -0+ and ip_ , Eq. (|2]) 
can be decomposed into the following coupled mode equa- 
tions: 
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J^+ + ^_) t (^++^_)](^ + + ^_)=0, (4) 

with l co h = koL co h quantifies the ratio between the co- 
herence length and the lattice constant. 

In passing, we remark that the above coupled mode 
equations are similar to the nonlinear coupled mode 
equations arising in the study of nonlinear wave prop- 
agation in one-dimensional periodic structures [l6[ . 

Equation (2) is equivalent to the Gross-Pitaevskii 
equation with an optical lattice potential. While trans- 
ort problems in BEC have been studied extensively 
the present study differs from the BEC systems in 
several ways. One of the important differences is that 
the non-equilibrium driving conditions arises much more 
naturally in the present system. Preparing a coherent 
source of BEC is much more difficult compared to a co- 
herent source of photons, i.e., a laser. Another is that we 
study transport properties in a sinusoidal lattice poten- 
tial equivalent to an optical lattice potential, whereas in 
the BECs, transport studies are often limited to a single 
or double well potentials, or for a periodic delta function 
potential that are analytically tractable. 

It is perhaps helpful here to mention that, unlike other 
nonlinear systems mentioned earlier, two periodicities are 
present in the current optical system. The first one is due 
to counter-propagating control fields creating stationary 
excitations [23, [28[ and the second one, created by the 
modulated atomic density, provides the periodic poten- 
tial in the NLSE shown in Eq. (pQ). 

Following Hafezi et. al. [10|, we specify the boundary 
conditions ip+(z = 0) = a and ip-(z = d) = 0, where 
a is the field strength for the driving laser impinging 
from the left and d = koL is the dimensionless length 
of the waveguide. This represents a situation depicted 
in Fig. [TJa) where the only input field impinges on the 
waveguide from the left. The transmission spectra, i.e., 
the transmittivity as a function of the probe field detun- 
ing e, is investigated for both linear and nonlinear regimes 
in the semiclassical limit and especially the effects of the 
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lattice potential is explored. Detuning the probe field 
corresponds to detuning the quantum fields which can 
be easily implemented in experiments. Consequences of 
this detuning is discussed in Appendix iBl 



III. LINEAR REGIME 

If the nonlinearity is assumed to be zero, Eq. (j2j) is 
equivalent to Mathieu's differential equation with well 
known solutions. In this section, we plot the transmission 
spectra obtained from these solutions. The lattice poten- 
tial introduces bandgaps and bound states and shifts the 
resonance peaks in the transmission spectrum. Inciden- 
tally, locations of these band gaps can be calculated eas- 
ily from the imaginary part of the Mathieu Characteristic 
Exponent as shown in Appendix O 

Before we show you the transmission spectrum, we 
would like to point out that the linear case can be 
achieved without the 4th level in Fig. [T|d, which makes 
experiments easier [28|. In the linear case, one is able 
to study experimentally a classic transmission problem 
of a quantum particle through a potential. To put it an- 
other way, using the periodic optical medium it is possi- 
ble to build a tunable optical experiment to simulate the 
Schrodinger equation in a controlled environment. The 
shape of the potential can also be varied to test other 
classic text-book problems by changing the atomic den- 
sity distribution. 





Figure 2: Transmission spectra for l co h — 0.25, d — 3tt, f3 — 
12tt/50 « 0.75, V = 0(a) and V = 1.0 (b). The bound state 
and the first band gap is clearly shown in (b). 

Figure [2] shows the transmittivity as a function of 
the input field detuning for two different values of V 
for l co h = 0.25 with total loss determined by f3 = 
(d/l co h)(T/\Ao\ = 127r/50 « 0.75. Bare optical parame- 
ters corresponding to these effective parameters are given 
in Appendix O where it is shown that the magnitude of 
V can be changed by tuning A p . The values were care- 
fully chosen such that the phenomena depicted in this 
work are within reach of experimental investigation. 

When V = 0, there are perfect transmissions at res- 
onances eo = (nir) 2 for n = 1,2,3,... and at an addi- 
tional point e~o = 1. These transmission resonances are 
an artefact of an effective cavity produced by counter- 
propagating control fields. For non-zero V, a continuum 
exists after 2V (maximum potential height) as expected, 



but there are additional resonances for e < 2V signifying 
bound states that exist due to the potential. The energy 
of the first bound state can be readily calculated within 
the perturbation theory as V — V 2 /8, which matches the 
starting point of the first resonance in the figure. Fur- 
thermore, the location of other resonances can be approx- 
imated from a generalization of the free field resonance 
points d\ft = nir to d^/e — V ± V 2 /8 = nir, such that 
only the quasi-kinetic energy part is used to determine 
the resonance condition. This approximation works well 
for e ^> V, i.e., when the effect of the potential is only 
perturbative and the solution is expected to be close to 
the free field case. It can be easily seen that this expres- 
sion matches the resonances for e > 2V. 



IV. NON-LINEAR REGIME 

In this section, we study the effects of the nonlinearity 
and the lattice potential on the transport properties. As 
remarked earlier, the nonlinear transport problem with 
a potential has been studied in detail in various phys- 
ical systems. For example, Paul et. al. and Rapedius 
et. al. have studied similar nonlinear barrier transmission 
problems for BECs, where they showed the bistability in 
the transmission spectrum [5|-[Sl- However, the present 
system offers unique advantages over the previous works: 
tunability of the parameters, large nonlinearity and read- 
ily available coherent sources. There are many interest- 
ing phenomena in the non-linear regime, but in this work 
we focus on the bistability (or multistability) as a typical 
phenomenon in nonlinear systems. We choose the param- 
eters used in the previous section and Re (x) = ±0.02 and 
±0.05. The sign and the magnitude of x can be tuned by 
varying the single-photon detuning A p as shown in the 
Appendix [Al 
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Figure 3: Transmission spectra for V — 0.5, l co h = 0.25, a = 
0.5, d — 37T, j3 ~ 0.75 and three different values of |%|. (a) For 
repulsive nonlinearity (x — 0.02), one can see the shift of the 
transmission peaks to the right (black curve). For attractive 
nonlinearity (% = —0.02), the shift is in the opposite direction 
(red curve), (b) Effect of the nonlinearity on transmission 
spectra is shown for x — 0-05 (black curve) and % = —0.05 
(red curve). The blue curve depicts the transmittivity for 

x = o. 

Figure [3] shows the effects of repulsive and attractive 
nonlinearity on the transmission spectrum. For repulsive 



nonlinearity (black curves), one can see the shift of the 
transmission resonances to the right, whereas for attrac- 
tive nonlinearity (red curves) the shift is in the opposite 
direction. The nonlinear shift induces strong bistability 
in the first peak. The suppression of the transmission 
peaks is due to the presence of nonlinear loss absent in 
the linear case. Increasing the magnitude of nonlinear- 
ity of course causes larger nonlinear shift and more loss. 
In Fig. [3fb) the shift is sufficient to make the second 
resonance of the repulsive case overlap with the third 
resonance of the attractive case. Further increase in the 
nonlinearity for the attractive case results in the trans- 
mission peaks disappearing to the left. This phenomenon 
signifies the shift of the transmission peak from the res- 
onance state to a bound state at e < [H . 

Fixing x an d increasing the input field strength a has 
similar effects on the transmission spectrum to fixing a 
and increasing \- This is expected because the effec- 
tive nonlinearity determining the nonlinear shift is de- 
termined by xl^l 2 , where the over bar denotes averaging, 
and |^ 2 of course depends on a. 

Having observed the nonlinear shifts and bistability, 
we investigate the effect of the lattice potential on these 
nonlinear effects. We concentrate on the repulsive case 
for this purpose. Figure |4] depicts the scenario for three 
different values of V = 0.1, 0.5 and 1.0. 




Figure 4: Transmission spectra for l co h = 0.25, a = 1.0, 
X — 0.02, d — 3tt and fi ~ 0.75. (a), (b) and (c) correspond 
to V = 0.1,0.5 and 1.0 respectively, (c) shows an interesting 
behaviour near e = 2. A horizontal peak, near the bandgap 
region, arises due to an interplay between lattice potential 
and nonlinearity. 

As V increases, the transmission spectrum shifts to- 
wards the higher value of e. The constant shift is due to 
the constant part of the potential which shifts the energy 
and the reduced transmittivity of the third peak results 
from the development of a band gap. This shift and 
the development of the band means that the resonance 
peaks get pushed towards each other as V increases. Near 



V = 1.0, the second and the third peak merges and a pe- 
culiar feature is observed between e = 1.5 and e = 2.0: a 
horizontal peak near the bandgap region appears, as seen 
in Fig. lU^c). As far as the authors are aware, this hori- 
zontal peak has not been observed before and is a unique 
feature arising due to an interplay between the nonlin- 
earity and the band-gap effect of the lattice potential. 
We have verified that this feature is more pronounced 
for extreme parameters and its development is under in- 
vestigation. 

Finally, we remark that not all of the multiple solu- 
tions are stable. Normally the solution corresponding 
to the negative slope is unstable and the system shows 
bistable behaviour. To show which part of the transmit- 
tivity curve is occupied one has to study the dynamics of 
the system as in [5j for example. A hysteresis behaviour 
can then be observed as a relevant parameter is varied 
in time. In this way a connection to the swallow-tail 
behaviour discussed in the BEC literature can be made 
through a general argument given by Mueller [29| . 



V. CONCLUSION 

We have studied a nonlinear optical system based on 
EIT and showed that this system is useful for studying 
linear and nonlinear transport properties. In the linear 
case, the system provides a clean and tunable labora- 
tory realization of a textbook quantum transport prob- 
lem. In the nonlinear case, the system allows the study 
of the transport problem in an optical lattice, similar 
to that studied in the context of a BEC in an optical 
lattice. Transport is a lot easier to investigate in the 
proposed optical system, due to the ease of preparing a 
coherent source compared to the BEC case. Transmis- 
sion spectrum have been studied for various regimes of 
parameters and properties such as band gap and multi- 
stable behaviour have been demonstrated. The authors 
believe the prevalence of the nonlinear Schrodinger equa- 
tion studied in this work and the non-equilibrium driven- 
dissipative nature of the system makes the work interest- 
ing to a broad audience. 
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Appendix A: Experimental feasibility 

Various effective parameters appearing in Eq. (|2]) and 
below can be written in terms of bare optical parameters 



5 



as 
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(Al) 
(A2) 
(A3) 
(A4) 

(A5) 



where A = Vt 2 /(Vt 2 — 5Ao/2), Y\d denotes the sponta- 
neous emission rate into the guided modes and V the 
total spontaneous emission rate from the excited states 
| b) and \d) assumed to be the same for simplicity. The 
detunings are shown in Fig. [I] (b) and the last line fol- 
lows immediately from l co h = k®/ (ArriRVg). Thus, to get 
a large nonlinearity l co h has to be small, because to sup- 
press losses it has been assumed that |A P | ^> T > Y\d- 

In the main text, we fix the value of the nonlinearity 
X- Once this is done, the value of V can be changed 
according to the ratio: 



V 

X 



2Q 2 



n\ 1 
k A' 



(A6) 



We will assume the parameters mostly similar to those 
used in all the frequencies are in units of the spon- 
taneous emission rate which we take to be T = lOMhz 
for simplicity. S = -0.01, A = -50, T 1D = 0.2, 
n 10 7 m _1 , ni lO 6 ™ -1 , k lO 4 ™ -1 . With these 
values we obtain 



lcoh — 0.25, 



V 



1 



10 (ft 2 -0.25)' 

1 \ 0.2A W A 2 



2A P J A 2 + 1/4' 

Values of the real part of x (Xr) an( ^ V as functions of 
A p and Q are plotted in Fig. [5] 



Appendix B: Probe field detuning 

With the numbers used in the previous section, we de- 
duce that Er w 5 x 10 _3 r, assuming that v g = 100ms -1 . 
Because the recoil energy is several orders of magnitudes 
smaller than the one photon detunings, one can safely ig- 
nore the contributions to the one photon detunings due to 
the probe field detuning. One cannot, however, arbitrar- 
ily increase the two-photon detuning to explore the spec- 
trum, because if the two-photon detuning is too large, the 




Figure 5: Plots of (a) the real part of the interaction param- 
eter XRi as a function of one photon detuning A p . (b) The 
lattice depth as a function of Rabi frequency Q. Here, both 
A» and Q are in terms of F. 



two-photon resonance condition breaks and EIT works no 
more. The two photon resonance condition is Q? ^> \ST\, 
where S is the two photon detuning. Thus, one will have 
to modify the frequency of the control laser to keep S 
within the EIT regime. We will assume that this is the 
case in this work. 



Appendix C: Mathieu characteristic exponent 

As mentioned in Sect.[TTI]the locations of the bandgaps 
can be found from the Mathieu characteristic exponent 
(MCE), where the imaginary part is non-zero at band 
gaps. This follows from the fact that all Mathieu func- 
tions F(a, q, z) can be written in the form e w P(a, q, z), 
with P{a 1 q 1 z) periodic in z. As an example we plot the 
transmittivity and the imaginary part of the MCE for 
V = 2, d = 3tt and l co h = 0.1 in Fig. [6] shows. Finite 
length of the potential causes the small non-zero trans- 
missions seen at the band gap. 




Figure 6: Transmittivity (blue solid line) and the imaginary 
part of Mathieu Characteristic Exponent (red dotted curve) 
for V — 2, d = 3tt and l co h =0.1. MCE predicts the band 
gaps quite well even for the relatively short waveguide length 
used. 
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